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Abstract. Formal asymptotic expansions have long been used to study the 
singularly perturbed Allcn-Cahn type equations and reaction-diffusion systems, 
including in particular the FitzHugh-Nagumo system. Despite their successful 
role, it has been largely unclear whether or not such expansions really represent 
the actual profile of solutions with rather general initial data. By combining 
our earlier result and known properties of eternal solutions of the Allcn-Cahn 
equation, we prove validity of the principal term of the formal expansions for 
a large class of solutions. 



1. Introduction 

In this paper, we study the behavior of solution u £ of an Allen-Cahn type equation 
of the form 

1 



(P e ) 



Ut = Au + -~(f(u) - sg £ (x, t, u)) in 17 x (0, oo) 

— = onfflx (0, oo) 

dv 



u(x, 0) = uq(x) in . 

and also that of a reaction-diffusion system of the form 



u t = Au ■ 



1 



(f(u) + sf 1 (u,v)+e'fl(u,v)) 



(RD £ ) I 



vt = DAv + h(u, v) 
du dv 
dv dv 
u{x, 0) = uq(x) 

v(x,0) = v (x) 



in il x (0, oo) 
in Six (0, oo) 
on dVl x (0, oo) 
in il 
in f2 , 
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where f(u),g £ (x,t,u),fi(u,v) and /f(it,«) satisfy the conditions to be specified 
later, and £ is a positive parameter. A typical example of (RD e ) is the FitzHugh- 
Nagumo system: 

u t = Au + \ (f{u) - efi(u) - ev) in Six (0, oo) 

v t = DAv + au — /3v in X (0, oo) 
(FHN e ) I du dv 

v ; ' — = — = on oil x (0,oo) 
av ov 

u(x, 0) = uq(x) in Q, 

v(x, 0) = vq(x) in Q . 

It is well-known that the solution u £ (x,t) of the above systems develops a steep 
transition layer, which converges to a "sharp interface" as e — > 0. To study such 
a sharp interface limit, formal asymptotic expansions of u £ are commonly used to 
discover, formally, the law of motion of the limit interface. Then, based on these 
expansions, one can construct sub- and super-solutions or use some approximation 
argument to prove the convergence of the transition layer — or the front — to the 
sharp interface, thereby establishing rigorously that the limit motion law agrees 
with what is anticipated from the formal asymptotics. 

However, this standard approach only tells us that the transition layer of u £ is 
confined within a relatively narrow zone — of thickness o(l) or sometimes even 
0(e) — around the limit interface, but it does not say much about whether or not 
the actual transition layer really possesses a robust profile that matches the formal 
asymptotics. Known answers to this question are mainly concerned with solutions 
whose initial data already has a well-developed transition layer. The case of more 
general solutions has largely been unexplored. Our goal is to provide an affirmative 
answer in this direction: we shall prove that, for e sufficiently small, the solution 
u £ — with rather general initial data — of both (P e ) and (RD e ) possesses a profile 
that agrees with the principal term of the formal expansion. 

1.1. Notation and assumptions. The notation and assumptions stated below 
strictly follow those in [2]. 

In problems (P E ) and (RD e ) above, is a smooth bounded domain in M. N (N > 2) 
and v denotes the outward unit normal vector to dfl. The nonlinearity is given by 
f(u) := —W'(u), where W(u) is a double-well potential with equal well-depth, 
taking its global minimum value at u = a± . More precisely we assume that / is C 2 
and has exactly three zeros a_ < a < a + such that 

(1.1) /'(a±)<0, /'(a) > (bistable nonlinearity) , 
and that 

(1.2) / f(u)du — (balanced case). 

J a_ 

In the Allen-Cahn equation (P e ) we allow the balance of the two stable zeros a_ 
and a+ to be slightly broken by the function —sg £ (x, t, u) defined on fl x [0, oo) x K. 
We assume that g £ is C 2 in x and C 1 in t, u. and that, for any T > there exist 
•6 G (0, 1) and C > such that, for all (x, t,u) G O x [0, T] x E, 

e\A x g £ (x,t,u)\ + e\g £ t (x, t, u)\ + \g £ u (x,t,u)\ < C, 

\\g £ (-,-,u)\\ ... i+* _ < C. 
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Moreover, we assume that there exists a function g(x,t,u) and a constant, which 
we denote again by C, such that 

(1.3) \g e (x,t,u) - g(x,t,u)\ < Ce , 

for all small e > 0. In [2], we also assumed = on dfl x [0, 00) x K, though this 
last condition is only for technical simplicity. Note that these conditions, except 
the last one, are automatically satisfied if g e is smooth and independent of e. 

In the reaction-diffusion system (RD £ ) we assume that fi(u,v), f^(u,v) are C 2 
functions and that /f , along with its derivatives, remain bounded as e —> 0. We also 
assume that D > and that h(u, v) is a C 2 function such that, for any constants 
L, M > 0, there exists a constant Mi > M such that h(u, —Mi) > > h(u, Mi) for 
M < L. These conditions enable us to construct a family of invariant rectangles 
mentioned in Remark 1. 

In the FitzHugh-Nagumo system, we assume that fi <G C 2 (M.) and that a and /3 
are given positive constants so that (FHN e ) becomes a special case of the reaction- 
diffusion system (RD e ). 

To complete the picture we need to specify conditions on the initial data. We 
assume that uq and vo belong to C 2 (f2). We define the "initial interface" Tq by 

(1.4) r := {x e : U (x) = a} , 

and assume that To is a C 3+ ^ hypersurfacc (0 < 1} < 1) without boundary such 
that 

r CC fl and Vu (x) • n(x, 0) 7^ if x G T , 

uq > a in Sid , uq < a in J7q , 

where denotes the region enclosed by Tq and Qq the region enclosed between 
dfl and To, and n(x,0) denotes the outward unit normal vector at x € Tq = dflQ . 
Let us emphasize that we do not assume that the initial data uq of if already has 
well-developed transition layers depending on e, in which case the validity of the 
formal expansions is more or less known (see subsection 1.2 for more details). 

Remark 1 (Time-global smooth solutions). Under the above assumptions, it is 
classical that (P e ) has a uniformly bounded smooth solution u s that exists for all 
t > 0. As for (RD e ), the same can be shown for e > small enough, by using the 
method of invariant rectangles (see e.g. [2] for details). 

1.2. Known results for the singular limit. We present here a brief overview of 
known results. Hcuristically, in the very early stage, the diffusion term is negligible 
compared with the reaction term. Hence, in view of the profile of /, the value of u E 
quickly becomes close to either a+ or a_ in most part of f2, creating a steep interface 
(transition layers) between the regions {u e w a_} and {u e ~ a + } (Generation of 
interface). Once the balance between diffusion and reaction near the transition 
layers is established, the interface starts to propagate in a much slower time scale 
(Motion of interface). The interface obeys a certain law of motion, which is to be 
investigated. 

A first step to understand this motion is to use (inner and outer) formal asymp- 
totic expansions of u £ . This was performed in the pioneering work of Allen and 
Cahn [3] and, slightly later, in Kawasaki and Ohta [18], who revealed that the in- 
terface motion involves curvature effects. Using such arguments, one discovers that 
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the sharp interface limit of (P e ) obeys the following law of motion: 
V n = —(N — 1)k + c / g(x,t,r)dr on T t 

J a- 



(P°) 



r, 



t=0 



r . 



where T t denotes the limit sharp interface at time t > 0, V n is the normal velocity of 
r t in the exterior direction, k the mean curvature at each point of r t , cq a constant 
determined straightforwardly from / via 



CO 



V2 / (W(s) - W(a-)) 1/2 ds 



where W(s) := — J s f(r)dr. As long as the solution T t of (P°) exists, we denote by 
il^ the region enclosed by T t , and by fi^ the region enclosed between 5£1 and T t - 
Also we define a step function u(x, t) by 



(1.5) 



u(x, t) 



in f2 t 
in fl+ 



to which u £ is formally supposed to converge as e — > 0. As regards (RD e ), the limit 
problem is found to be 



(RD°) < 



V n = -(JV-I)k-co / fi(r,v)dr on f, 



St = DAv + h(u, v) 
r *lt=o = r ° 

*L = 

ov 

^ ?)(x,0) = u (a;) 



in f2 x (0, oo) 

on 3f2 x (0, oo) 
in n , 



where u is the step function defined in (1.5). This is a system consisting of an 
equation of surface motion and a parabolic partial differential equation. Since u is 
determined straightforwardly from r t , in what follows, by a solution of (RD°) we 
mean a pair (r,{>). 

Remark 2 (Local smooth solutions for the limit problems). Under our assumptions, 
there exists T max > such that (P°), respectively (RD ), possesses a unique smooth 
solution T = Uo<t<T m »«= (rt x {i}), rcsp. (r, v) = (Uo<«T"'« (r t x {t}), v). For more 
details we refer to [2] and the references therein, in particular [13], [14], [12]. In the 
sequel we select any < T < T max and work on [0, T]. 

Numerous efforts have been made to rigorously prove the convergence of (P e ) 
and (RD e ) to (P°) and (RD ), respectively. Concerning the Allen-Cahn equation, 
let us mention the work of de Mottoni and Schatzman [20] (generation of interface 
via sub- and super-solutions) and [21] (motion of interface via construction of and 
linearization around an ansatz) or that of Bronsard and Kohn [10] (motion of inter- 
face via T-convergence) . Chen [11, 12] has established an 0(e\ lne|) error estimate 
between the location of the actual transition later and the limit interface, both for 
scalar equations and systems for rather general initial data. More recently, in [2], 
the present authors improved this estimate to 0(e). More precisely, they show that 
the solution u £ develops a steep transition layer within the time scale of 0(e 2 \ lne|), 



VALIDITY OF EXPANSIONS IN ALLEN-CAHN AND FITZHUGH-NAGUMO 



5 



and that the layer obeys the law of motion that coincides with the formal asymp- 
totic limit (P°) or (RD°) within an error margin of O(e). Let us also mention that 
there are results of much finer error estimates in the literature (see for instance [8]), 
but those results are concerned with very specific initial data which already have 
nice transition layers consistent with formal asymptotics (hence dependent on e). 

As mentioned before, in most of the aforementioned works, approximate solu- 
tions or sub- and super-solutions are constructed by roughly following the formal 
expansions. Our goal is to investigate the actual validity of such expansions for 
solutions u e with rather general initial data. 

Remark 3 (Viscosity framework). Since the limit problem may develop singu- 
larities in finite time, the classical framework does not always allow to study the 
singular limit procedure for all t > 0. Nevertheless — as far as the Allen-Cahn 
equation is concerned — following [17], [15] one can define a limit problem for all 
t > that generalizes (P°) in the framework of viscosity solutions. In this setting 
we refer to [16] (convergence of Allen-Cahn equation with prepared initial data to 
generalized motion by mean curvature), [4], [6] (generalizations), [22, 23], [7], [5] 
(not well-prepared initial data), [1] (fine convergence rate). 

2. Main results 

We start by giving an outline of the formal asymptotic expansions mentioned 
before. See [2, Section 2] for more details. Let u e be the solution of (P e ), and 
r = Uo<t<T(rt x {t}) be the solution of the limit geometric motion problem (P°). 
We define the signed distance function to T by 



(2.1) d(x,t) 



— dist(a;, T t ) for x £ £l t 
dist(£, r t ) for x G Clf . 

Then, near T, we make a formal inner expansion of the form 
(2.2) u e (x, t) = U (x, t, ) + eUi (x, t, < ^^- 



(2.3) 



Some normalization conditions and matching conditions (with the outer expansion) 
are also imposed. By plugging the expansion (2.2) into (P e ), we discover that 
Uo(x,t, z) = Uo(z), where Uo(z) is the unique solution (whose existence is guaran- 
teed by the integral condition (1.2)) of the stationary problem 

U " + f(U a ) = 

Uq(— oo) = a_ , Uq(0) = a , Uq(oo) = a + . 

This solution represents the first approximation of the profile of a transition layer 
around the interface observed in the stretched coordinates. Next the solvability 
condition for the equation involving U\ provides the law of motion (P°) for the 
limit interface T, which, in turn, determines the term U\. 
It is then natural to wonder if the ansatz 

u^x,t) = U (^)+eU 1 (x,t, d{x > t) 



e J \ e 

is really a good approximation of the profile of the solution u e . Note that the 
convergence results mentioned in subsection 1.2 do not answer this question; indeed 
those results simply show that the level surface of the solution u £ 

(2.4) := {x G fl : u e {x,t) = a] 



() 
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converges to the sharp interface (T t )o<t<T, which is a solution of (P°), and that 
(2.5) lim u £ (x,t) -- 



for x G ri t 
for x <E f2^~ , 



without clarifying the validity of (2.2). Our main result Theorem 2.1 below provides 
a first answer in this direction. In the sequel we define 

(2.6) t^/VrVllnel, 

which is the time needed for the transition layer of u e to become fully well-developed 
(see Lemma 3.1). We define the signed distance function associated with T e by 



(2.7) d E {x,t):= 



-dist(x, rf ) if u € (x, t) < a 
dist(x,r|) if u £ (x,t) > a. 



Note that this definition of d £ is consistent with that of d in (2.1) in view of (2.5). 

Theorem 2.1 (Validity for Allen-Cahn). Let the assumptions of subsection 1.1 
hold (in particular the initial condition uq is rather generic). Let u e be the smooth 
solution of Allen-Cahn equation (P e ). Fix fx > 1. Then the following hold. 

(i) If e > is small enough then, for any t € [fit e ,T], the level set T E is a 
smooth hypersurface and can be expressed as a graph over T t . 

(ii) 



(2.8) lim sup 



fj,tc<t<T,xen 



U e (x, t) — Uq 



d e {x,t) 



0. 



where d e denotes the signed distance function associated with T £ 
(iii) There exists a family of functions 

(2.9) 6 E : U <t<T(r t x {*}) ->R (0 < e « 1) 

whose L°° -norms remain bounded as e — ¥ 0, such that 

'd(x,t) - e6 e (p(x,t),t) 



(2.10) lim sup 



ntz<t<T,xen 



U S (x, t) — Uq 



= 0. 



where d denotes the signed distance function associated with T and p{x, t) 
denotes a point on T t such that dist(x,T t ) = \\x — p(x,t)\\. 

Note that p(x,t) is an orthogonal projection of the point x onto r t , which is 
uniquely defined in a small tubular neighborhood of T t since T t is a smooth solution 
of (P°). Note also that the presence of the perturbations —£0 e (p(x,t),t) cannot be 
avoided since it reflects the small difference between d(x,t) and d £ (x,t). 

Let us mention that the validity of higher order terms of the formal expansions 
(for generic solutions) is still unknown. 

Our next theorem provides similar estimates for the reaction-diffusion systems. 

Theorem 2.2 (Validity for the reaction-diffusion system). Let the assumptions of 
subsection 1.1 hold, and let (u 6 ,v 6 ) be the smooth solution of the reaction- diffusion 
system (RD e ). Fix fi > 1. Then the same conclusions as in Theorem 2.1 hold, 
with d being the signed distance function associated with (T,v), which is the smooth 
solution o/(RD°) on [0,T]. 
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3. Proof of the main results 

The proof of Theorem 2.1 relies on the following two results: 

(a) the level set T £ is approximated by the interface T t by order 
0{e) ([2], see subsection 3.1 of the present paper) , 

(b) any eternal solution that lies between two planar waves is actu- 
ally a planar wave ([9], see subsection 3.2 of the present paper), 

combined with a rescaling argument. 

3.1. Thickness of the layers: the refined 0(e) estimate. We quote a result 
which is valid for both (P E ) and (RD £ ). 

Lemma 3.1 ([2, Theorem 1.3 and Theorem 1.11]). Let rj be an arbitrary constant 
satisfying < T) < min(a — — a). Then there exist positive constants Sq and 

Co such that, for all e £ (0, £o) and for all /'(a) _1 e 2 | lne| = t E < t < T, we have 

(3.1) \u s (x,t)-a±\<v if xent\Af Co e(r t ), 

where J\f r (Tt) ■ = {x £ £1 : dist(x,T t ) < r} denotes the r -neighborhood ofT t . This 
implies in particular that T E C Afc oe (^t) for all t e < t < T , hence 

(3.2) \d £ (x, t) — d(x, t)\ < CqS for all (x,t) £Q x [t £ ,T] , 0<e «1. 

3.2. Eternal solutions and planar waves. We recall that a solution of an evo- 
lution equation is called eternal (or an entire solution) if it is defined for all positive 
and negative time. We follow this terminology to refer to a solution w{z, r) of 

(3.3) w T = A z w + f(w) , zeR N ,reR. 

Stationary solutions and travelling waves are examples of eternal solutions. Cru- 
cial to our analysis is a recent result of Berestycki and Hamel [9] asserting that "any 
planar-like eternal solution is actually a planar wave" . More precisely, the following 
holds (for zeM^we write z = (zW, • • • , z^)). 

Lemma 3.2 ([9, Theorem 3.1]). Let w(z, r) be an eternal solution of (3.3) satisfying 

(3.4) liminf inf w(z, r)>a, limsup sup w(z, r)<a. 



where z' := (z^\ • • • , z^ N 1 '). Then there exists a constant z* £ R such that 

w(z,t) = U (z^ - z*), zeR N lT eR. 

3.3. Proof of (ii) in Theorem 2.1. In what follows we fix \i > 1 and an arbitrary 
constant T\ with T < T\ < T max (see Remark 2). Obviously the conclusion of 
Lemma 3.1 remains valid if T is replaced by T±. Assume by contradiction that (2.8) 
does not hold. Then there is 77 > and sequences Ek I 0, tk £ [fit £k ,T], x^ £ ft 
(k = 1, 2, ...) such that 

'd^{x k ,t k y 



(3.5) 



u £k (x k ,t k ) - U Q 



> 2?y. 



In view of (3.1)-(3.2) and C/o(±oo) = a±, for (3.5) to hold it is necessary to have 

(3.6) d(x k ,t k ) = 0(e k ) , asfc^oo. 

If u 6k (x k ,t k ) = a, then this would mean that x k £ rf^, in which case the 
left-hand side of (3.5) would be (since Uq(0) = a), which is impossible. Hence 
u Ek (x k ,t k ) a. By extracting a subsequence if necessary, we may assume without 
loss of generality that u £k (x k ,t k ) — a has a constant sign for k = 0, 1, 2, Since the 
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sign of this quantity is irrelevant in the later argument, in what follows we assume 
that 

(3.7) u £k (x k7 t k )>a (k = 0,1,2,...), 
which then implies that 

d e *{x k ,t k )>0 (k = 0,1,2,...). 

Since the sequence (x k ) remains close to T tk by (3.6), and since r tfc is uniformly 
smooth for k = 0, 1, 2, . . ., each x k has a unique orthogonal projection p(x k , t k ) £ 
r tfc . Let y k be a point on T e t * that has the smallest distance from x k . If such a 
point is not unique, we choose one such point arbitrarily. Then we have 

(3.8) u E «(y k ,t k ) = a (k = 0, 1, 2, . . .), 

(3.9) d £ *(x k ,t k ) = \\x k -y k \\, 

(3.10) u e *(x,t k )>a if \\x-x k \\ < \\yk-xk\\ , 

x k - p k _l r tfc at p fe e r tfc , 

where p/j := p(x kl t k ). Furthermore, (3.6) and (3.2) imply 

(3.11) ||a*-Pk||=0(eiO, \\Vk-Pk\\ = 0(e k ) (fc = 0,1,2,...). 
We now rescale the solution u e around (p/b,ifc) and define 

(3.12) w k (z, t) := u £k (p k + e k TZ k z, t k + e 2 k T) , 

where TZ k is a matrix in SO(N, R) that rotates the z'"' axis onto the normal at 
Pk € r ffc , that is 

ft fe : (0,...,0,l) T ^n(p fc ,t fc ), 

where ( ) T denotes a transposed vector and n(p,t) the outward normal unit vector 
at p G T t . Since T t (hence the points p k ) is uniformly separated from <9f2 by some 
positive distance, there exists c > such that w k is defined (at least) on the box 

B k :=\(z,T)eR N xR : ||z|| < -, — (/j — l)/'(a) _1 | lns^l < r < -^—5— 

I £ k £ k 

Since u e satisfies (P e ), we see that w k satisfies 

(3.13) w k = A z w k + f{w k )-e k g e x{p k + e k K k z,t k + e 2 k T,w k ) in B k . 
Moreover, if (z,r) 6 -B fc then t Efc < t k + e\t < T\. Therefore (3.1) implies 

j d(p k +e k K k z,t k + e 2 k r) < ~C e k =S> w k (z, r) < a_ + r\ , 
I d(p fc + s k K k z,t k + e 2 k r) > C e k => w k (z,r) > a + - rj , 

so long as (z, t) 6 _B fc . Since the rotation by lZ k of the z^ N > axis is normal to r tfc at 
Pk '■= p(x k ,t k ), and since the curvature of T t is uniformly bounded for < t < T, 
we see from (3.14) that there exists a constant C > such that 

(3.15) z W < -C=> w k (z,T) < a_ + r), Z W >C^w k {z,T)>a+-r 1 , 

so long as (z,r) e B fc and ||zj| < e^ 1 . 

Now, since w fc solves (3.13), the uniform (w.r.t. k > 0) boundedness of w k and 
standard parabolic estimates, along with the derivative bounds on g £ , imply that 
w k is uniformly bounded in C,^ 7 ' +2 (-B 1 ). We can therefore extract from (iu fc ) a 
subsequence that converges to some w in Cf^B 1 ). By repeating this on all B k , 



(3.14) 
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we can find a subsequence of (w k ) that converges to some w in C ; 2 Q '* (R N x R) (note 
that U k > a B k = R N x K). Passing to the limit in (3.13) yields 

w T = A z w + f(w) on R N x K . 

Hence we have constructed an eternal solution w(z, r) which — in view of 

(3.15) — satisfies (3.4). Lemma 3.2 then implies that 

(3.16) w(z,t) = U (z^ -z*) 

for some z* £ E. 

Now we define sequences of points (zfe), (zit) by 

z fe := — TZ7 1 (x k -Pk) , h ■= —R-^iVk -Pk) ■ 

By (3.11), these sequences are bounded, so we may assume without loss of generality 
that they converge: 

Zk Zoo , z k -t Zoo , as k V OO . 
By the definition of the z coordinates, Zoo must lie on the z^ N ' axis, that is, 

z oo = (0,...,0,*£°) T . 
It follows from (3.8) and (3.10) that 

(3.17) w(z oo ,0)=a, w(z,0)>a if \\z — Zoo\\ < ||5oo - 2<»|| • 

Note that by (3.16), the level set w(z, 0) = a coincides with the hyperplane z^' = 
z*, and recall that Uq' > 0. Therefore, in view of (3.16) and (3.17), we have cither 
Zoo = Zoo, or that the ball of radius \\zoo — Zoo|| centered at Zoo is tangential to the 
hyperplane z^ N ' = z* at Zoo- This implies that Zoo, as well as Zoo, must also lie on 
the z( N ) axis. Therefore 

z oo = (0,...,0,z*) T , 
and the inequality w(zoo,0) > a implies that Z6o^ > z*. On the other hand (3.9) 

implies d 6k (x k ,t k )/e k = \\x k - yk\\/£k = \\zk - z k \\ ->• ll^oo - Zoo|| = 4^' - z*. The 
assumption (3.5) then yields 

= \\w(zoo,0)-U (z^ ~z*)\\ 

= || lim u e "(x k ,t k )-U ( lim £^5*1**)) || 

A;— >-oo \ / 

> 277 . 

This contradiction proves statement (ii) of Theorem 2.10. 

3.4. Proof of (i) and (iii) in Theorem 2.1. The proof of (i) below uses an 
argument similar to the proof of Corollary 4.8 in [19]. Fix fi > 1. For a given 
77 £ (0,min(a — a_,a+ — a)) define £0 > and Co > as in Lemma 3.1. Then we 
claim that 

(3.18) liminf inf \7u £ (x, t) ■ n(p(x, t), t) > , 

V e^O xeMc oe (T t )^t'<t<T 

where n(p, i) denotes the outward unit normal vector at p £ T t . Indeed, assume 
by contradiction that there exist sequences e k 4 0, t k £ [pt £k ,T], x k £ J^fc e k {^t k ) 
(k = 1, 2, ...) such that 

Vu £k (x k ,t k ) ■ n{pk,tk) < 0, 



10 



M. ALFARO AND H. MATANO 



where pt = p{xk,tk). By rescaling around (jPk,tk) and using arguments similar to 
those in the proof of (ii), one can find a point with \z£ | < Co such that 

which contradicts the fact that Uq' > and establishes (3.18). Since, in view of 
Lemma 3.1, Tf C Afc s(^t), the estimate (3.18) implies that Vu e (x,t) ^ for all 
x E Tf; hence by the implicit function theorem, is a smooth hypersurface in a 
neighborhood of any point on it. The fact that can be expressed as a graph over 
r t also follows from (3.18). This proves statement (i) of Theorem 2.1. 

Finally, statement (iii) follows immediately from statements (i), (ii) and (3.2). 
This completes the proof of Theorem 2.1. 

3.5. Proof of Theorem 2.2. As shown in [2, Section 7], the behavior of u e in 
the system (RD e ) can be treated as a special case of (P e ), by regarding v e as a 
given function and using a contraction mapping theorem. Thus the conclusion of 
Theorem 2.2 follows directly from Theorem 2.1. 
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